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Abstract 

The quantum renormalization group method is applied to study the quantum criticality and 
entanglement entropy of the ground state of the Ising chain in the presence of antisymmetric 
anisotropic couplings and alternating exchange interactions. The quantum phase transitions can 
be characterized by the discontinuity in the second derivative of the energy of renormalized ground 
state. The phase diagram is obtained by the critical boundary line. The first derivative of entangle- 
ment entropy also diverges at the same critical points after enough iterations of the renormalization 
of coupling constants. The antisymmetric anisotropy and alternating interaction can enhance the 
renormalized entanglement via the creation of quantum fluctuations. The scaling behavior of the 
derivative of the entropy around the critical points manifest the logarithm dependence on the size 
of the spin system. 
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I. INTRODUCTION 



The resource of entanglement can play an important role in the quantum information 
processing such as quantum teleportation and algorithms for quantum computation jl, 2]. 
The entanglement is quantum correlation which can describe the property of the ground state 
for various spin systems 3J. It has been justified that there exist the close relations between 
the entanglement and the quantum phase transitions in condensed matter physics j^-lO]. 
The quantum critical phenomena happen at zero temperature when the spin couplings and 



external fields are varied in the vicinity of critical points [11] . As is known, the quantum 
phase transitions can be represented by the discontinuities in the derivatives of some order 
parameters. However, for no a priori knowledge of suitable order parameters, some special 
physical quantities can be used to universally identify the quantum criticality at the ground 



state 



12j,|l3|. The ground-state energy is considered as one characterizing the quantum phase 



transitions for some spin systems 



14 



151 ] . So far some measures of quantum entanglement 



have been found out to be relevant to the quantum criticality for many spin systems 



16 



181 ] . During the theoretical investigation of low-dimensional spin systems, the numerical 



methods of density-matrix renormalization group and exact diagonalization have received 



much attention 



19 



211 ] . In many materials, the magnetic properties were studied by the 



spin models with the antisymmetric anisotropic interaction such as Dzyaloshinskii-Moriya 
interaction 



22 



23] . Recently, the quantum renormalization group approach hasprovided a 



convenient way to analyze the quantum phase transitions for such systems 24M27T] . By this 
method, the renormalized ground states can be described in terms of matrix product states 
The quantum phase diagram can be approximately obtained by the enough iterations 



of the renormalization of couplings and external fields. Meanwhile, rich quantum criticality 
can also be embodied by the bond-alternating spin models 29( . These respects motivate us 
to investigate the effects of antisymmetric anisotropy and alternating spin interactions on 
the entanglement and quantum phase transitions. 

The paper is organized as follows. In section II, the quantum renormalization group 
method is introduced and used for the Ising chain with antisymmetric anisotropy and alter- 
nating exchange interactions. The renormalized ground state can be obtained analytically. 
In section III, the quantum phase transitions are approximately characterized by the en- 
ergy of the renormalized ground state. The boundary line of critical points determines two 
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different phases. The effects of antisymmetric anisotropic interactions and alternating ones 
on the entanglement entropy are investigated. The divergence of the derivative of entropy 
occurs at the same critical points. The finite size scaling behavior of the entanglement are 
also studied. A simple discussion concludes the paper. 



II. RENORMALIZATION OF THE SPIN MODEL 



The implementation of quantum renormalization group is the perturbation method to 



30|, 



reduce the degrees of freedom of the spin model. In the Kadanoff's block approach 
the original Hamiltonian of the model can be divided into the block Hamiltonian and the 
interacting Hamiltonian. The block Hamiltonian is the sum of some commuting items which 
act on different blocks. The effective Hamiltonian can be obtained by mapping the original 
Hamiltonian to the low energy subspace via the projector consisting of some lowest energy 
eigenstates of the block Hamiltonian. 

The Hamiltonian of the Ising chain with the antisymmetric anisotropic and alternating 
interactions can be written by 

JV 

H=jJ2ll~ (-iy\]S*S* +1 + D ■ (Si x S i+1 ). (1) 

i=l 

The periodic boundary condition of iV+ 1 = 1 is considered. Here Si = (Sf, Sf , Sf) denotes 
the Pauli operator acting on the z-th spin and Sf\ t (4)) = ±|| t (4))- The case of J > rep- 
resents the nearest-neighboring antiferromagnetic coupling and < A < 1 describes the rela- 
tive strength of alternating exchange coupling. The antisymmetric anisotropic interaction is 
described by the vector of the Dzyaloshinskii-Moriya interaction, D = (D x , D y , D z ), arising 
from the spin-orbit coupling. In the following analysis, the case of D x = D y = 0, D z = D is 
taken into account. 

According to the renormalization group method, the original Hamiltonian is divided into 
two parts 

L 

H = H° + H I = J2h° k + hi k+1 . (2) 
fc=i 

The part of H° = X/A=i^fe ^ s ^ e block Hamiltonian where L = N/3 represents the 
number of blocks and each item has three sites, = JYm=iI^ + (~^) k+l ^\^ik^i+ik + 
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D(Sf k Sf +lk — Sf k Sf +lk ). Here S" k (a = x,y,z) denotes the a-component pauli oper- 
ator on the Z-th spin of the fc-th block. The other part of H 1 represents the inter- 
acting Hamiltonian between two neighboring blocks, H 1 = J2k=i k+i where h k k+1 = 
J{[1 - (-l) k \]Sl k Sl k+1 + D(Sl k Sl k+l - SfAVi)}" We need the degenerate lowest en- 
ergy eigenstates \ip k ) of h k to construct the projector = | lUXV'fc I + I •U'fcXV'fc I where 
| f|~ (JJ-)) is the renamed basis in the effective Hilbert space for each block. When k is odd, 
the analytic expression of \ip k ) can be given in the product Hilbert space of a block 

i^ + ) = «im)+*i m)+c\ itt) 

\K) = a\m-ib\m)+c\m)- (3) 

Here the real parameters of a, b, c satisfy the condition of the normalization where a = 
bD/(X - 2e ),b = + 2D 2 (4e 2 + \ 2 )/(4e 2 - A 2 ) 2 . Here e is the scale of the lowest 

energy by the parameter J and calculated by the minimal root of the equation 8cq + 4€q — 
2(A 2 + 2D 2 )e - A 2 = 0. For the simplest example of A = 0, e = -(1 + Vl + 8L> 2 )/4. 
When k is even, the degenerate lowest energy eigenstates have the similar expression in Eq. 
(3) and the real parameters b even = b, a even = — c, c even = —a. The effective Hamiltonian 
mapped onto the low energy subspace can be expressed by 

L L 

He ff = ([{P k )H([[Pi). (4) 

k=l k=l 

By the effective spin operators fk for the k-th. block, the definite form in the effective Hilbert 
space can be written by H ef f = J' J2k=A^ ~ (~^) k M T k T k+i + P > '( T k T t+i ~ T k T k+i) where the 
operator r z = |(| ft)(fT | — | JJ-)(JJ- |). For the case of k = odd, the relations between the spin 
operator S^k = (Sf k , S y k , Sf k ) and the effective one fk = (r k , r k , r k ) are given by 

PkS^kPl = {2abr y k , -2abr x k , (1 - 2c 2 )t*) 
P k S 2 ,kPl = (2acr y k , 2acr x k , (1 - 2b 2 )r z k ) 

PkS^Pl = {2bcr y k , -2bcT%, (1 - 2a 2 K) . (5) 

The renormalization of the couplings satisfy the equation 

/^(l-2 -)(l-2^ g ^- (1 _ 2 4 ;; ( f_ 2e2) . (6) 

In general, after the n-th renormalization group iteration step, the chain of 3 n+1 sites can be 
represented by the effective three-site Hamiltonian with the renormalized couplings. Accord- 



ing to the renormalization of ground state given by the expression of Eq. (3), the quantum 
phase transitions and the entanglement properties can be studied. 



III. PHASE DIAGRAM AND ENTANGLEMENT ENTROPY 

The quantum phase transitions at the ground state are driven by the quantum fluctua- 
tions. It is the key to select the suitable order parameters. To obtain the phase diagram, we 
choose the ground-state energy as one universal measure. From the renormalization equa- 
tion of couplings in Eq. (7), it is seen that the quantum phase transitions will happen by the 
change of the antisymmetric anisotropic interaction D and relative strength of alternating 
spin exchange A. For a certain value of A, the obvious discontinuity in the second derivative 
of ground-state energy can be shown at the quantum critical point D c in Fig. 1(a). It 
means the feature of the second-order quantum phase transition. With the change of the 
alternating interaction A, the phase diagram is obtained in Fig. 1(b). It is shown that two 
different phases are divided by the critical points D c ~ y/1 — A 2 . 

It is interesting to investigate the entanglement properties of the spin system in the two 
different phases. Using the renormalization group method, we calculate the entanglement 
between some degrees of freedom and the rest of the system. The renormalization of entan- 
glement in the large enough iteration step can be used to estimate the large size behavior 
of the spin system. As one successful measure of entanglement, the block entropy can be 
implemented to describe the global property of the entanglement. The density matrix of the 
renormalized ground state is p — IV'oXV'oI where |^ ) — Wt) or IV'o") ^ s expressed by Eq. 
(3). Therefore, the reduced density matrix for the middle site can be obtained by tracing 
over the rest sites 



The entropy E for the middle site can represent the entanglement between the degrees of 
freedom of the middle site and ones of the rest sites. The expression of E is given by 



In the first step of renormalization group, the chain of 3 2 sites can be described by the 
effective Hamiltonian i? e //. The change of entanglement entropy E is depicted as a function 
of the antisymmetric anisotropic coupling D and alternating exchange interaction A in Fig. 




(?) 



E = -(l-b 2 )\og 2 (l-b 2 )-b 2 \og 2 (b 2 ). 



(8) 
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2. It is shown that the values of E are increased with D and A. The Dzyaloshinskii-Moriya 
interaction in the XY plane can induce the planar quantum fluctuations which enhance the 
quantum correlation of the ground state. The existence of alternating couplings can spoil the 
form of the Neel ordered state which has no entanglement. The values of the entanglement 
tend to one maximum when A is close to one. 

To analyze the connection of the entanglement to the quantum phase transition, we need 
study the scaling property. By means of increasing the iteration step, the entropy E can be 
calculated for the large-size spin system. From Fig. 3, it is seen that the rate of the change of 
entanglement entropy will be increased as the size becomes large through the iteration step. 
In the phase / for the case of D < D c , the entropy is always less than one. The maximal 
entanglement entropy of E = 1 occurs in the phase II for the case of D > D c . After the large 
enough iteration step, the appearance of nonanalytic behavior of the entanglement entropy 
verifies the feature of the quantum phase transition. In this spin system, the measure of 
entanglement entropy can be applied to estimate the quantum criticality. 

Furthermore, it is also necessary to evaluate the singularity of the derivative of the renor- 
malized entanglement entropy. As depicted in Fig. 4(a), the drop of the first derivative 
of entropy will become more pronounced with the increase of the spin sites. The minimal 
values of — |^ at some points D = D m are decreased rapidly with the size. The points of D m 
for a large iteration step are closely near to the critical ones D c . The scaling of the points 
D m is also shown in Fig. 4(b). The logarithm dependence on the sites N is numerically 
given by D m = D c — iV~ 7 where the parameter 7 = 7(A) is increased with the decrease of 
the alternating coupling A. Through the numerical fit, 7(1) ~ 0.46 and 7(0.5) ~ 0.52. This 
respect manifests that the divergence of the derivative of the entropy appears more rapidly 
for the smaller values of alternating interaction A in the condition of the same size. 

IV. DISCUSSION 

The entanglement property and quantum phase transitions in the Ising chain with the 
Dzyaloshinskii-Moriya interaction and alternating coupling are studied by the method of 
quantum renormalization group. Based on the analytic renormalization equation, the renor- 
malized ground state can be obtained. The ground-state energy is considered as one universal 
measure of quantum phase transitions. The nonanalytic behavior of the second derivative 
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of the renormalized ground-state energy can justify the feature of second-order phase tran- 
sition. The critical boundary line divides the ground state into two different phases. It 
is demonstrated that the antisymmetric anisotropy and alternating couplings can help for 
the increase of the entanglement entropy. The discontinuities in the entanglement entropy 
happen at the critical points which are same to those provided via the order parameter of 
the ground-state energy. With the decrease of the alternating coupling, the divergence of 
the derivative of the entropy appears more rapidly in the vicinity of the quantum criticality. 
In this one-dimensional spin system, the renormalized entanglement entropy is the efficient 
measure of the quantum phase transition. 
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Figure caption 

Figure 1 

(a). The second derivative of the energy of the renormalized ground state diverges at 
the critical point D c ~ 0.85 in the iteration step n = 9 when A = 0.5. (b). The dash line 
denotes the critical boundary characterizing two different phases / and II. 

Figure 2 

The renormalization of the entanglement entropy E is plotted as a function of the anti- 
symmetric anisotropy D and alternating coupling A in the iteration step n = 1. 
Figure 3 

The evolution of the renormalized entanglement entropy as the size is shown for A = 0.5. 
The square data denote the case of the step n = 1, the circle ones represent the case of 
n = 5 and the triangles are the data for the case of n = 9. 

Figure 4 

(a). The negative values of the first derivative of the entropy are calculated when A = 0.5. 
The dot line describes the case of iteration step n — 1, the dash line denotes the result of 
n = 5 and the solid one represent the case of n — 9. (b). The data characterize the finite 
size scaling property of the singular points D m where the values of — |^ are minimal. The 
squares represent the case of A = 0.5 and the circles denote the case of A = 1. 
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